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Abstract. Heavy quark production in high-energy proton-nucleus (pA) collisions is 
described in the framework of the Color Glass Condensate. k± factorization is broken 
even at leading order albeit a more general factorization in pA holds at this order in 
terms of 2, 3 and 4 point correlators of Wilson lines in the nuclear target. The x- 
evolution of these correlators is computed in the large A and large N mean field limit 
of the Balitsky-Kovchegov equation. We show results for heavy quark production at 
RHIC and LHC energies. 



1. Introduction 

Heavy flavor production in high energy processes is usually described in the collinear or 
k± factorization frameworks of perturbative QCD. Reactions with nuclear targets allow 
one to study new important features of multiple scatterings and parton saturation, which 
is characterized by the scale Qj^Jl]. For a large nucleus at higher energy, this scale can 
be large and comparable to the heavy quark mass tuq, so that Aq CD -C m 2 Q < Q 2 sA <C s. 
In this situation, realized at RHIC and definitely at future LHC, we can investigate 
saturation physics through heavy quark production with the generalized framework of 
the Color Glass Condensate (CGC) [2, 3, 4, 5, 6, 7]. 

In nucleus-nucleus collisions, heavy quarks are produced by initial hard interactions 
and subsequently propagate in the hot flowing medium which modifies the pj_ spectrum 
of the quarks and the quarkonium yield[8]. In order to quantify these in-medium effects 
in hot matter, proton/deuteron-nucleus (pA) collisions play a decisive role providing the 
baseline of the initial nuclear effects. Here, we will compute heavy quark production 
in high energy pA collisions using the CGC framework that incorporates both multiple 
scattering and small x evolution effects. 



2. Quark pair cross section in the large N limit 



In pA collisions, which is a prototype of a dilute-dense system, the particle production 
cross-sections are known analytically to the leading order in the strong coupling constant 
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a s and the proton source p p , but to all orders in the dense nuclear source g 2 pA — 0(1) [2]. 
The cross-section of the quark pair production is written in the CGC framework as 



d(7 HQ S 



d 2 p ± d 2 q ± dy p dy q (2tt) 6 C f 

X [ ^±^lM <f>f9 2 (k 2± ,k ± ) Vpm (k 1± ) , (1) 
k 2 ±,k± 

where ki± = p ± + q ± — k 2 ±, and fci )2 j_ and yi )2 = ln(l/xi i2 ) are the momenta and 
rapidities of the gluons coming from the proton and the nucleus. A shorthand notation 
S(fcij_, fc 2 ±, k±) stands for the matrix element squared, whose explicit form can be found 
in Refs. [2, 3, 4]. In the leading twist approximation on both the proton and nuclear 
sides, this S reduces to the LO hard matrix element in the k± factorization formalism[9]. 
The (p pm (ki±) is the fc^-dependent gluon distribution in the proton, which is related to 
the usual gluon distribution function via: 

i r Q2 

— J dk 2 ± <p P ,y(k ± ) = xG p (x, Q 2 ) . (2) 

The multi-gluon correlations in the large nucleus are encoded in the expression (1) as 
the 3-point correlator (large N) 

I 2 N 

4>™(l±, k ± ) = nR 2 A ^- S Y (k ± ) S Y (l ± - k ± ), (3) 

where Sy{k±) is the Fourier transform of the forward scattering amplitude of the right- 
moving dipole S Y (x ± ) = (l/N)tr(U(x ± )W (0)) Y with 

f +oo 



U(x±) = "P+exp 



/+oo 1 
dz + -2 P a A (z + ,x ± ) t a 
oo V I 



(4) 



Here t a is the SU(iV) generator in the fundamental representation. 

On the proton side, the collinear approximation is more appropriate because the 
typical k± should be small - O(Aqcd) as long as x\ is not too small, especially in the 
forward rapidity (large xi) region. As noted in Refs. [9, 2] the collinear limit on the 
proton side is well-defined thanks to a Ward identity: one obtains (see Fig. 1 (a)) 

dUqq a 2 



d? Pl _d 2 q ± dy p dy q 4(2n)*C F 

x / ^J^lM ^(fc 2± ,fc ± ) x 1 G p (x u Q 2 ), (5) 



where fc 2 _i_ = Pj_ + Q±, and 

-i(u u \- v [ d6i E (k 1± ,k 2 ±,k ± ) 

^ " J = ^oj ^ jr ' < 6) 

In this work, we numerically solve the Balitsky-Kovchegov (BK) equation for 
Sy(k±) [10, 11] with the McLerran-Venugopalan (MV) model[12] initial condition 
at xo=0.01, and compute the 3-point correlator 0^ 2 (fc 2 j_, k±) (Fig. 1 (b)). This 
prescription is justified in the large N and large A limit. For large-x extrapolation, 
see Ref. [4]. 
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(b) *■ < a u ' Q s 2 = 2.0 GeV 2 




Figure 1. (a) Collinear limit on the proton side. Blob stands for the interactions with 
the nucleus including multiple scatterings, (b) Nuclear gluon distribution obtained 
from the BK equation with the MV model initial condition. 
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Figure 2. (a) Open charm spectrum at y/s = 200 GeV at y=0, 2 calculated with 
a s =0.2 and Q^ Q= 2GeV 2 . (b) Nuclear modification factor of Q 2 sA0 =2 GeV 2 to 0.22 
GeV 2 for D meson at -^5=200 and 8500 GeV. In the latter case, results of the BK 
evolution with a s =0.1 and 0.2 are compared. 



3. Open heavy flavor 

In Fig. 2 (a), we show the single charm spectrum from eq. (5) with Q 2 sA =2 GeV 2 at 
x = 0.01 and the BK evolution with fixed a s =0.2. CTEQ6LO[13] is used for the proton. 
The D° spectrum is obtained after convolution with Peterson's fragmentation function 
D(z) (e = 0.05 and z = pn±/Pc±)- Preliminary data are taken from STAR [14, 15]. 
We see that after the convolution the D° spectrum is a little harder at lower p± than 
the data while it fits in the higher p± region. Note that the b-quark contribution is 
not included in our calculation. Fig. 2 (b) presents the nuclear modification factor of 
Q 2 sA (x = 0.01)= 2 to 0.22 GeV 2 , at y/s= 200 GeV (y=0, 2) and 8500 GeV (y=0, 2, 
4). We see at RHIC energy the Cronin peak moving to higher p± with the rapidity y 
because the charm pair production is dominated by the moderate x 2 region. Instead, 
at LHC energy we see larger suppression with y; Note that fixed coupling a s =0.2 yields 
too large A of Q 2 s (x) ~ (l/x) A compared with HERA analysis. The saturation effects 
can be seen even in the B meson productions at this energy 
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4. Quarkonium 

The J ftp production cross-section is computed for Q 2 sA0 =2 and 0.22 GeV 2 in the color 
evaporation model in Fig. 3, where we assume that the bound state formation takes place 
outside of the nucleus, and is thus not modified compared to pp. The J/t/> absorption in 
the target is not necessary. The cross-section is supressed (enhanced) in the small (large) 
p± region at RHIC energy due to the multiple scattering of the pair (Fig. 3 (a)) [16]. 
At forward rapidities dN/dy is more suppressed due to the multiple scattering and 
saturation, which may be compared with the data[17]. Note that in the hatched region 
y < 0.5 in Fig. 3 (b), dN/dy is largely determined by the larger x 2 part of 4>a,'L, which 
is extrapolated from the MV model. In Fig. 3 (c), we show the suppressions of dN/dy 
for the J/tjj and the charm quark at LHC energy. We expect that the difference between 
them represents the effect of the multiple scattering on the bound state formation. It 
is interesting to note that this difference depends on the rapidity y only weakly. 
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Figure 3. (a) J/ip spectrum for Q 2 sA0 =2 and 0.22 GeV 2 at ^fs~= 200 GeV. (b) 
Nuclear modification factor of and charm quark (a s =0.1 and 0.2). Hatched region 
correspoinds to x-i > 0.01. (c) The same as (b) at s/s= 8500 GeV. 
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